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ABSTRACT 

We present the first luminous, spatially resolved binary quasar that clearly inhabits an ongoing 
galaxy merger. SDSS J125455.09+084653.9 and SDSS J125454.87+084652.1 (SDSS J1254+0846 here- 
after) are two luminous z=0.44 radio quiet quasars, with a radial velocity difference of just 215km/s, 
separated on the sky by 21kpc in a disturbed host galaxy merger showing obvious tidal tails. The 
pair was targeted as part of a complete sample of binary quasar candidates with small transverse 
separations drawn from SDSS DR6 photometry. We present follow-up optical imaging which shows 
broad, symmetrical tidal arm features spanning some 75kpc at the quasars' redshift. Previously, the 
triggering of two quasars during a merger had only been hypothesized but our observations provide 
strong evidence of such an event. SDSS J1254-I-0846, as a face-on, pre-coalescence merger hosting 
two luminous quasars separated by a few dozen kpc, provides a unique opportunity to probe quasar 
activity in an ongoing gas-rich merger. Numerical modeling suggests that the system consists of 
two massive disk galaxies prograde to their mutual orbit, caught during the first passage of an active 
merger. This demonstrates rapid black hole growth during the early stages of a merger between galax- 
ies with pre-existing bulges. Neither of the two luminous nuclei show significant instrinsic absorption 
by gas or dust in our optical or X-ray observations, illustrating that not all merging quasars will be in 
an obscured, ultraluminous phase. We find that the Eddington ratio for the fainter component B is 
rather normal, while for the A component L/L^dd is quite (>3it) high compared to quasars of similar 
luminosity and redshift, possibly evidence for strong merger-triggered accretion. More such mergers 
should be identifiable at higher redshifts using binary quasars as tracers. 

Subject headings: black hole physics; galaxies: active; galaxies: interactions; galaxies: nuclei quasars: 
emission lines 



1. INTRODUCTION 

The origin, growth, and evolution of massive galaxies, 
and the supermassive black holes (SMBHs) that they 
host, represents a prime field of study in modern as- 
trophysics. We n ow know that galaxies re gularly in- 
teract and merge (jToomre k Toomrd flQTl) . and' that 
SMBH r eside in the centers o f most, if not all galax- 
ies (e.g., iRichstone et al.lll998[) . These two facts alone 
suggest that binary SMBHs should be commonplace. 
Of course, one or both of the SMBHs in a binary will 
only be detectable as quasars when they are actively 
accreting. One of the leading proposed mechanisms 
to trigger stro ng accretion (quasar) activity is galaxy 
mergers (e.g., iHernguistl Il989t iKauffmann fc Haehnelti 
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120001 : I Hopkins et al.l I2008L and references therein) , so 
merging galaxies with bin ary quasars should also be 
common. (Bcgclm an et al.1 [i980) first discussed binary 
SMBH evolution, from galaxy merger to coalescence, 
as an explanation for the form and motion of radio 
jets in Active Galactic Nuclei (AGN) The 'Final Par- 
sec Problem' (iMilosavljevic fc MerrittI l2003f l - whether 
the coalescence of a b i nary SMBH ultimately stalls 
(|Milosavlievic fc Me rritt 2001) , proceeds to rapid coa- 
lescence fe.g.. Escala et al. 2004), or in stead recoils or is 
ejected (e.g.. IMadau fc Ouataertil2004l ) - has important 
implications for the detection of gravitational waves and 
for the spin and demography of SMBHs. 

In a broader context, astronomers hypothesize that 
"feedback" - whereby dynamical interactions between 
galaxies trigger accretion onto their SMBHs - medi- 
ates the tight correlation between galaxy central black 
hole masses and t he velocity dispersions cr^, of galaxy 
bulge s {Mbh — CT-^ : lFerrarese et al.ll2000l : iGebhardt et al.l 
[2000 ). The resulting quasars grow in the galaxy cores 
until they blow out the very galactic gas that feeds 
them (e.g.. iGranato et al.ll2004[ ). choking off star forma- 
tion , and eventual l y leading to passive elliptica l galax- 
ies (jHoEMnilerS] ISqoI iKomend^i^ This 
feedback paradigm dovetails with cosmological models 
of hierarchical structure formation if quasar activity is 
induc ed by massive mergers (e.g., Wvithe fc Locb 200^ 
[20051) . Major mergers (i.e., those with mass ratios above 
^ 0.3; iSheiJ l2009f l between gas-rich galaxies most ef- 
ficiently channel large quantities of gas inwards, fos- 
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tering starbursts and feeding rapid black hole growth. 
Deep high resolutio n imaging of quasa r host galaxies 
(iBahcall et alj Il997t iGuvon et all l2006t iBennert et all 
I2008D shows strong evidence for fine structure and 
tidal tails expected from past gravitational interactions. 
Radio-quiet quasar hosts tend to be found in gas-rich 
galaxy mergers that form intermediate-mass galaxies, 
while radio-loud QSOs reside in massive early-type galax- 
ies, most of w hich also show signs of recent mergers 
or interactions (jWolf fc: She ini^ '2008*) . The far- infrared 
(FIR) emission of QSOs appears to follow a merger- 
driven evolution from FIR-bright to FIR-faint QSOs 
(|Veilleux et al.ll2009aD . 

The measured excess of quasars with < 40 kpc separa- 
tions fe.g.. lHennawi et al.ll2nnfil : IMvers et alllBoTl [2008h 
over the extrapolated large-scale quasar correlation func- 
tion may indeed be due to mutual triggering, but is also 
debated to ari se naturally from th eir locally overdense 
environments (jHopkins et al.|[2008f) . The dynamics and 
timescales of major mergers are therefore of the utmost 
interest. 

To date, the merger hypothesis is supported by find- 
ings of spatially-resolved binary active galactic nu- 
clei (AGN) in just a handful of z < 0.1 galax- 
ies with one or both of th e nuclei heavi l y ob- 
scured in X-rays (NGC62 40, iKomossa et"all 120031: 
Arp299, IZezas et al.l 120031: Mrk463. IBianchi et at] 
[2008), b y the unusua,l BL Lac- type obje ct OJ 287 
([Sillanpaa et al.l 119881 iValtonen et all 120091 ) . and per- 
haps by X- shape d m o rphology in radio galax ies (e.g. , 
IMerritt fc E kcrs 2Q0l lLhil[200l ICheundl2007l) . In ad- 
dition, COSMOS J100043. 15-^020637.2 is known to con- 
tain two AGN resolved at 0.5" (~1.8kpc) separation in 
HST/ACS imaging, which have a radial velocity differ- 
ence of Aw = 150 km/s , and appear to be ho sted by a 
galaxy with a tidal tail (jComerford et al]|2009f) . 

Interest in spatially unresolved systems (spectro- 
scopic binary AGN candidates) has surged of late, 
spawned largely from the troves of SDSS spectroscopy. 
The unusual sy s tem SPSS J153636. 22-1-1044127.0 
(|Boroson fc Laueil I2009f ) has a spectrum with two 
broad line systems separated b y At; = 3500 km/s , 
and also has a nearby radio (Wrobel fc Laoii l2009f) 
and optical ([Decarli et al.l 12009 ) counterpart. The 
physical natu re of this system has been heavily de- 
bated (e.g., [ Chornock ct al. 2010"; ILauer fc BorosonI 
120091: iTangfc Grindlav 2009 ). largely because spatially 
unresolved quasars with double-peaked broad emis- 
sion lines are a q uite common phenomenon (e.g., 
iStrateva et al.l I2003D . Another spatially unresolved 
system, SDSS J092712.65+294344.0, shows 2 broad 
and 1 narrow emission line system in its spectrum, 
with Av — 2650 km/s, sparkin g discussion abou t 
whether it is a chance super position (iShields et al.ll2009D . 
hosts a recoiling SMBH (jKomossa et al.l I2008D . or IS 
a bound binary SMBH inside a single narrow line re- 
gion (Bogdanovic ct al. 2009; Dotti ct al. 2009). Quasars 
with double-peaked narrow emission lines are relatively 
comm on (e.g., there ar e at least 167 such systems in the 
SDSS; ILiu et al.ll2009l). From the D EEP2 Galaxy Red- 
shift Survev. iComerford et al.l ()2009[ ) found more than a 
third of Type 2 AGN show [O HI] line velocities signif- 
icantly (50-300 km s~^) offset from the redshifts of the 



host galaxies' stars, arguing that the most likely expla- 
nati on is inspiral ling SMBHs in merger-remnant galax- 
ies. iSmith et all (|2009) find that about 1% of (21,592) 
quasars in SDSS DR7 have detectable double-peaked 
[O HI] emission line profiles. Only 2 of those appear to 
he spatially resolved, but the (single-fiber spectroscopy) 
sample selection is strongly biased against nuclei with 
separations greater than about an arcsecond. 

Spatially unresolved systems such as these are rela- 
tively easy to find in large spectroscopic samples. How- 
ever, because of the lack of spatial information, the veloc- 
ity offsets are open to a variety of interpretations depend- 
ing on the relative strength and velocity of narrow and/or 
broad emission line systems: small-scale gas kinematics, 
asymmetric or thermally inhomogeneous accretion disks, 
AGN outflows or jets, recoiling or orbiting SMBHs, or 
distu rbed or rotating narrow line regions (|Smith et al.l 
|2009( ). Furthermore, spectroscopic samples are biased 
against binary AGN that are very close in redshift (i.e., 
unresolved in velocity space), or those that are more 
widely separated on the sky. For instance, in the SDSS 
spectroscopic survey, the fiber diameter is 3", and the 
minimum separation of fibers on a plate is 55" on the 
sky. So, except in rare cases with multiple overlapping 
spectroscopic plates, any binary quasar with separation 
in between these two values could only be found from 
dedicated follow-up spectroscopy. 

Why are spatially-resolved active nuclei in mergers so 
rare? First, they may be heavily shrouded and there- 
fore only detectable as ultraluminous infrared galaxies. 
'ULIRGs' have bolometric luminosities rivalling quasars, 
and by some {HST I-band) estimates , as ra any as 40% 
retain double active nuclei (jCui et al.l 1200 It) . A binary 
fraction in local ULIRGs of at least 40% is also consis- 
tent with R- and K-hand ground-based data obtained 
under <, < 1" seeing and later confirmed with ff-band 
HST data (iVeilleux et al.l I2002L [200^ . Among dust- 



reddened quasars. lUrrutia et al 



( 20081 ) found that 85% 



show evidence of merging in images of their host galax- 
ies. Second, detectable mergers may be rare simply be- 
cause the lifetime of the resolvabl e-but-unmerged inter- 
acting phase i s extr emely short (jMortlock et al.l 119991 : 
iForeman et"all I2OO90 . Third, gas-rich major mergers 
should trace quasars, and therefore should mainly have 
occurred near t he "quasar epoch" a t higher redshifts 
(z ^1.5; e.g., IKhochfar fc Burkertl l200lt IWolf et al.l 
l2003t [Silverman et al.ll2005l ). where detection of extended 
host galaxy light is challenging. 

Th e prevailing vi ew in the li te rature (e.g.. IDiorgovskil 
19911: Kochanek ct~aL[ 119991: IMortlock et al.l 119991: 



Mvers et al.l 120071) is that the excess of quasars with 



small (< 40 kpc) separations is evidence for nuclear 
triggering in galaxies during dissipative mergers. Ac- 
cording to Hop kins et al. (2007), the excess measured 
clust ering (e.g., IHennawi et "all 120061: IMvers et al.] 120071 
I2008f ) indeed represents compelling evidence for the 
merger-driven origin of quasars. However, they also 
note that attaching all quasars to moderately rich dark 
matter environments in which mergers are most likely to 
occur is sufficient to explain the observed excess of bi- 
nary quasars at < 40kpc, even if they are not triggering 
each other in a bound orbit. That is, they just happen 
to be neighbors where the typical observed velocity 
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differences could represent ~ Mpc separations along 
the line of sight rather than dynamical velocities, and 
their properties should be statistically indistinguishable 
from those of single quasars. The discovery of binary 
quasars whose hosts are clearly interacting thus presents 
rare opportunities to study what merging/triggering 
really looks like, and allows for derivation of important 
quantities associated with the interaction. 

In the next section we discuss the discov- 

ery of SDSS J1254+0846, a pair of luminous quasars 
with nearly identical redshift, hosted by a galaxy 
merger. Unique among known spatially resolved sys- 
tems, SDSS J1254-f 0846 is at a scale which suggests it is 
an ongoing merger rather than a relaxed or remnant sys- 
tem. SDSS J1254-I-0846 can thus be used to help study 
boundary conditions for gas-rich galaxy mergers. The 
description and results of our deep follow-up multiwave- 
length imaging and spectroscopy follow in [JSj including 
estimates of the black hole masses and Eddington ratios. 
To verify the plausibility of the merger scenario and de- 
termine if it's possible to infer any properties of the in- 
teraction, we compare numerical iV-body simulations to 
the observed properties of the system in 5|4l In SJSJ we 
consider the hypothesis that the pair might be lensed, 
which we find to be extremely unlikely. We present our 
conclusions in [JBl Throughout, we assume the following 
cosmological parameters for distance-dependent quanti- 
ties: flm = 0.3, r^A = 0.7, and Hq = 72 km s'^ Mpc"\ 
which yields an angular size scale of 5.5 kpc/arcsec at the 
redshift of our system (z ^ 0.44). 

2. DISCOVERY OF SDSS J1254+0846 

The objects SDSS J125455.09-t-084653.9 
(SDSS J1254-I-0846 A hereafter) and 
SDSS J125454.87-I-084652.1 (SDSS J1254-I-0846 B 
hereafter) were targeted as a pair 

(SDSS J1254-H0846 hereafter) as part of a complete 
sample of binary quasar candidates with small trans- 
verse separations drawn from SDSS DR6 photometry 
(Myers et al. 2010, in preparation). A preliminary 
targeted follow-up campaign of such objects (for DR4) 
is discussed in Myers et al. (2008). Quasar candidates 
were selected as having g < 20.85 and either the 
"UVX quasar" and/or "low-redshift quasar" Bayesian 
classification flags set in the catalog of Richards et al. 
(2009)0 These cuts ensure a high efficiency of quasar 
pairs in the targets and a reasonably homogeneous 
sample over redshifts of 0.4 < z < 2.4. Pairs of quasar 
candidates were then followed up spectroscopically if 
they had an angular separation of 3" to 6". To extend 
the completeness of the sample as a function of comoving 
separation, the sample was also extended to pairs with 
separations of 6" to ^ 7.7" if neither component had a 
known redshift at z > 1.2. 

Following an extensive observational campaign with 
the R-C Spectrograph on the Mayall 4-m at Kitt Peak 
National Observatory and the Double Spectrograph on 
the 200-inch Hale telescope at Palomar Observatory, 
the sample of target quasar pairs from which SDSS 
J1254-I-0846 was drawn is now complete (again, see My- 
ers et al. 2010, in preparation). SDSS J1254-F0846 it- 
self was discovered to be a binary quasar on the night 

* uvxts=l OR lowzts=l 



of February 11, 2008 and confirmed as such on Febru- 
ary 12, 2008 at Kitt Peak (Myers and Hennawi observ- 
ing). As was the case for all candidate binary quasars 
observed on that run, a 1.5" by 204" long-slit set at 
the position angle of SDSS J1254-(-0846 was used, al- 
lowing both components to be simultaneously observed. 
The KPC-lOA grating was used yielding a resolution 
of 5A and a wavelength coverage of ~ 3800-7800A. 
As the seeing was ^ 1.5" on the nights in question, 
and the component separation of SDSS J1254-f0846 is 
3.8", care was taken when reducing the data to pre- 
vent the component spectra from merging. The com- 
ponent spectra were reduced iteratively using xidl's low- 
redux package (Hennawi, Buries, Schlegel & Prochaska; 
http : //www . ucolick . org/~ xavier/LowRedux/) with 
the procedure guided by hand using a boxcar extraction 
to ensure no overlap of the spectra. 

At 1.4 GHz, the host and its quasars are unde- 
tected, with flux less than 2.5 mJy at a resolution of 
45" (250kpc] in th e NRAO VLA Sky Survey (NVSS; 
iCondon et al.l[T998^ and less than 1 mJy at a resolutio n 
of 5" (28 kpc) in the F IRST survey C White et al.l[T997l) 1. 
Following llvezic et al.l (|2002,) , combining the 1 m Jy up- 
per limit with the «-band magnitudes of the quasars, both 
components are radio quiet. 

As the sample of binary quasar candidates from which 
it was drawn is now complete, SDSS J1254-I-0846 should 
be unremarkable. On the other hand, it is the lowest 
redshift binary quasar currently known and it has an un- 
usually low color similarity statistic (see Hennawi et 
al. 2006, Myers et al. 2008) of 0.2, meaning that the 
colors of its two components are practically identical. 
Across all 5 bands, their SDSS (PSF) magnitudes dif- 
fer by 2.28±0.08 (flux ratio 8.27±0.61), identical within 
the errors. Based on the {g — i) vs. redshift (Green et 
al. 2009) for SDSS quasars, the optical colors of these 
quasars are marginally (~ 0.2mag) blue relative to the 
mean, but consistent with those expected at this red- 
shift. The striking feature of SDSS J1254+0846 was dis- 
covered as we imaged this pair and several others in our 
Chandra/NOAO joint program (P. Green, P.I.) to ob- 
served binary quasars and their environments. Images 
we obtained on the nights of 18 March 2009 at NOAO's 
Kitt Peak Observatory with the MOSAIC imager on the 
4-m Mayall telescope (Barkhouse and Myers observing) 
immediately revealed bright tidal tails emanating from 
the quasar pair. The uniqueness of the system led us to 
procure further deep imaging and spectroscopy at other 
facilities. 

3. OBSERVATIONS 
3.1. Optical Spectroscopy 

We obtained deeper spectroscopy of both quasars on 
22 May, 2009 simultaneously through a single 0.9" slit in 
seeing of ^0.4" using the IMACS spectrograph on the 
Baade-Magellan telescope at the Las Campanas obser- 
vatory in Chile,. We centered the slit on QSO B, at a 
position angle of 61 degrees to include QSO A. We used 
the f/2 camera mode with a 300 lines/mm grism, giving 
a wavelength range of -4000 ~ 9600A and dispersion 
of 1.34A/pixel. We combined four exposures of 1200s 
each, and flux-calibrated using the original SDSS spec- 
trum, except for > 9200A where we used the standard 
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star LTT 3864, rescaled to match the SDSS spectrum, in 
the overlap region 8900 - 9050A. 

Figure [1] overlays the two spectra, with the B compo- 
nent scaled up and shifted for clarity. Despite the factor 
of ^11 difference in flux normalization, the redshift, con- 
tinuum and broad line shapes are all remarkably similar. 
The continuum ratio varies from about 11 near 5000A to 
about 10 near 8200A (observed frame). The discrepancy 
of 20 - 30% between flux ratios in our IMACS spec- 
troscopy and the SDSS imaging 8.27(±0.61) could be due 
to variability of either QSO component between the two 
epochs. This is consistent with the somewhat larger dis- 
crepancy in the blue, since QSO variabilit y is known to 
incre ase towards shorter wavelengths (e.g.. lWilhite et all 
12005% Some of the difference could also be due to slit 
position and alignment. 

The most striking difference is in the equivalent widths 
of the narrow emission lines. All the forbidden lines of 
[NeV], [Oil], [OIII] are relatively much stronger in B 
(larger equivalent widths) . To test whether the spectrum 
we have extracted for QSO B is contaminated by scat- 
tered light from A, we extracted a spectrum on the other 
side of QSO A at a distance equal to the separation of A 
from B. No significant spectral features are detectable, 
so we conclude that contamination of the B spectrum by 
scattered light from A is negligible. 

Figures [5] and [3] highlight the regions around H/3 
and H(5, respectively. The A spectrum looks some- 
what smoother because of its higher signal-to-noise ratio 
(S/N). The residuals plotted in the lower panels, from 
simple scaled subtraction with no velocity shift, are di- 
rect evidence that the redshifts of the two quasars are 
virtually identical. Separate cross-correlations of the two 
spectra against the SDSS median cornposite quasar spec- 
trum using IRAF xcsao (jKurtz et al.lll992| ). and exclud- 
ing telluric line regions and CCD artifacts yield za = 
0.43919±0.00028(i?=10.d3) and zb = 0.440095±0.00011 
(i?=15.0). From direct cross-correlation of the A and B 
spectra we find a velocity offset for A-B of -215±100km/s 
{R—6.7), consistent with essentially no significant veloc- 
ity difference. 

The residuals in Figures [2] and [3] also illustrate that the 
most significant differences between the spectra are in 
the narrow line components. However, some differences 
are also evident in the broad line regions. The spectro- 
scopic differences effectively preclude an interpretation 
of the pair as possibly lensed (see § [5] below for further 
discussion). Although it is difficult to tell from plots of 
this rescaled format, the A/B luminosity ratio is much 
smaller between the two components in low-ionization 
forbidden lines (ratio 1.6±0.2 in the [Oil] emission line) 
than in their continuum emission. While the [O II] emis- 
sion cannot be assumed to be a pure indica tor of star 
formation rate in the presence of an AGN (|Yan et al.l 
l200l . a larger fraction of [Oil] emission should arise 
from star formation in the host galaxy, so the smaller 
A/B ratio is evidence that the two nuclei probably share 
a host. 

We note that if these nuclei were spatially unresolved 
but with the same velocity difference and flux ratios, the 
system would not be detected in a spectroscopic SDSS 

^ R is the ratio of the correlation peak to the amplitude of the 
asymmetric noise. 



search for binary quasars such as that of Smith et al. 
(2009). Direct addition of the two components' spectra 
(as if in a single aperture) results in a rather normal- 
looking quasar spectrum - the B component is easily 
subsumed in the A spectrum, merely highlighting the 
very peaks of the narrow lines. On the other hand, if 
the tiny velocity difference is due to the serendipitously 
small angle of our sightline along the orbital axis, then a 
different projection would increase the observed velocity 
difference, as would a bound system with smaller physi- 
cal separation between the nuclei. 

3.2. Optical Imaging 

On the same night 22 May, 2009 at the Baade-Magellan 
telescope, we obtained 20 minutes of imaging (4 expo- 
sures of 300sec) in both Sloan r and i band. The IMACS 
f/2 camera has 0.2" pixels, and the seeing was 0.4". We 
subtracted the CCD bias level and flattened the fleld 
response using averaged dome projector flat images in 
each fllter as usual. Images were then co-added with the 
SWarp package (Bertin 2006, ver. 2.17.6), and object 
detection and measurement were made with SExtractor 
(Bertin & Arnouts 1996). Photometric calibration was 
performed using dereddened ma gnit udes from SDSS DR7 
for matching objects in the fleldPI The r-band image in 
Figure H shows the 2 bright nuclei of SDSS J1254-f0846 
and two symmetrical tidal tails spanning some 75 kpc at 
the quasar redshift. 

3.3. Chandra X-ray Observations 

We obtained X-ray images of the quasar pair with 
the Chandra X-ray Observatory on February 23, 2009 
at the ACIS-S aimpoint for 16ksec. The X-ray com- 
ponents are well-resolved by Chandra, and correspond 
closely (< 0.2") to their SDSS counterparts. To avoid 
cross-contamination, we extracted the X-ray photons 
from apertures corresponding to 90% of the counts (for 
1.5 keV). The NE (SW) components yield 1869 (381) net 
counts, respectively, in the 0.5-8 keV range. We fit an 
X-ray power-law spectral model 

N{E) =AE-^ X cxpi-N^^'aiE) - Nlja{E{l + Zabs))] 

to the counts using the CIAO tool Sherpa, where 
A is the normalization in photons cm~^ sec~^ 
keV~^ and a(E) is the absorption cross-s ection 
([Morrison fc McCamm(ml[T983t IwTlms et al.l[2?)00l ). We 
fix N^'^^ at the appropriate Galactic neutral absorbing 
column 1.9x 10^°atoms/cm~^, and include an intrinsic 
absorber with neutral column Nfj at the source red- 
shift. We group counts to a minimum of 16 per bin and 
fit using the statistic with variance computed from 
the data. The best-fit model for both components is 
F — 2.0 (with 90% confidence uncertainties of 0.05 and 
0.2 for the NE/SW components, respectively), and only 
upper limits to any intrinsic absorption {N"^*^ < 2.7 and 
7.2 xlO^° atoms/cm^ ^, respec tively). These valu es are 
quite typical of SDSS quasars (|Green et al.l 120091 ). The 
2keV A/B fiux ratio is 4.9, somewhat less than the op- 
tical flux ratio. The X-ray to optical ratio is often pa- 

1" We compare SDSS model_Mag to SExtractor MAG.AUTO val- 
ues. 
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rameterized by the X-ray-to-optical spectral slope cto J^ . 
which is 1.41 for A and 1.37 for B. X-rays in quasars 
become weaker relative to optical emission as luminosity 
increas es, and both t hese quasars fall a,long t he expected 
trends (jSteffen et al.ll2006l:lGreen et "alll2009[) . Statistical 
tests have shown that the correlation is weaker with red- 
shift, so that the aox{L) relationship is not a secondary 
effect of quasar evolution combined with the strong L — z 
trends of flux-limited quasar samples. While anecdotal, 
SDSS J1254-t-0846 confirms in a single system that the 
observed aox trend with luminosity in quasars is followed 
even by quasars at the same epoch and in the same large- 
scale environment. Both members of the pair are slightly 
X-ray bright for their estimated UV (2500A) luminos- 
ity. Our images suggest that the merger is essentially 
face-on between massive disk galaxies that are close to 
co-planar. If the quasar accretion disks are reasonably 
well-aligned with their galactic host diskj*^. than our 
sightline may simply avoid the extinction and redenning 
associated with a large angle to the line of sight. 

3.4. Radio Observations 

We observed SDSSjl254-H0846 with the VLA near 
transit on UT 2009 September 25 and 27, with net ex- 
posure times of 2350 s and 7053 s respectively, using the 
DnC configuration under NRAO proposal code AG826. 
We chose center frequency 8.4601 GHz (8.5 GHz here- 
after) with a bandwidth of 100 MHz for each circular 
polarization. Observations were phase-referenced to the 
calibrator J1254+1141 whose positional accuracy was 
less than 2 mas. The switching angle was 3° and the 
switching time was 240 s. Observations of 3C 286 were 
used to set the amplitude scale to an accuracy of about 
3%. The data were calibrated using the 2009 December 
31 release of the NRAO AIPS software. Each day's visi- 
bility data for SDSS J1254+0846 were concatenated and 
the AIPS task imagr was used to form and deconvolve 
a Stokes / image. Natural weighting was used to opti- 
mize sensitivity, giving an angular resolution at FWHM 
of 9.4" times 7.2" elongated at position angle -29°. One 
source was detected and an elliptical-Gaussian fit to it, 
yielding the following integrated flux density, position, 
and 1-dimensional position error: S = 0.26 ± 0.03 mJy, 
a( J2000) = 12''54'"55!08, (5(J2000) = +08°46'53"9, and 
fVLA — C.'S. The flux density error is the quadratic 
sum of the 3% scale error and the fit residual. The po- 
sition error is the quadratic sum of a term due to the 
phase-calibrator position error (less than 0.002"), the 
phase-referencing strategies (estimated to be 0.1"), and 
the signal-to-noise ratio (S/N) (0.3"). The source was 
unresolved and, given the modest S/N data, the corre- 
sponding diameter is less than the geometric-mean beam 
width, 8.2". 

The 8.5 GHz emission from SDSS J1254-I-0846 has a 
radio luminosity, Lfj — vLi,, of 1.5 x 10^° ergs s~^ and is 
unresolved, with a diameter of less than 8.2" (45 kpc). 
This scale encompasses the inner portions of the host 
galaxy merger, plus quasars A and B. Higher-resolution 
imaging using the Expanded VLA (jPerlev et all |2009[) 
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In ra dio galaxies, there is some general evidence against such 
alignment l lSchmitt et al.ll2002l) . 



is needed to localize the emission from each quasar. In 
the interim, we note that our Chandra data on quasars 
A and B impli es a 0.2-20 keV lumi nosity Lx of 7.9 x 
10^4 ergs s-i. ILaor fc BehaH (12001 propose that both 
active stars and radio quiet AGN owe their radio emission 
to similar coronal processes, following Lr/Lx ^ 10^^, 
where Lr = vL^, at 5 G Hz. For a spectral slope of —0.5 
(|Kellermann et al]|1994[ ). the observed 8.5 GHz luminos- 
ity corresponds to 1.1 x lO''" ergs s~^ at 5 GHz. Thus 
Lr/Lx 1.4 X 10~^, just the ratio expected for the 
radio-quiet quasars A and B. This simple, testable ar- 
gument suggests that the 8.5 GHz emission arises from 
both quasars A and B, without substantial contribution 
from the extended host galaxy. 

3.5. Black Hole Mass and Eddington Ratio 

Given our high quality spectra, we can estimate the 
black hole mass and Eddington ratios for each quasar. 
Most AGN black h ole mass estimators derive from rever- 
beration mapping ijPetersonI [19931 : iWandel et all [19991 . 
whereby time delays t between continuum and broad 
emission-line variations are used to deduce the size of 
the broad line region (BLR). For single-epoch optical 
spectra, the continuum luminosity La(5100A) can be 
used as a surrogate for the broad line region rad ius R 
(|Koratkar fc Gaskelll [l99l l : IKaspi et al.l[2000ll2005l) . By 
combining r with the emission-line width (most directly 
using H/3), a virial mas s for the black hole can be esti- 
mated (e.g., Vestcrgaard fc Pete rson 2006). 

We use the splot task within IRAF to deblend narrow 
and broad H/3 above the continuum in each quasar spec- 
trum, and we correct our H/3 line width measurements 
for spectral resolution and narrow line contamination fol- 
lowing Peterson et al. (2004) . For quasar A, we measure 
a H/3 FWHM of 2904±200 km s'^ . The log of the contin- 
uum luminosity at 5100A is 45.422 in erg/s. If we assume 
the bolometric correction o f 9.26 for 5100A luminosity 
from [Richards et all (|200g). then logLBoi = 46.39. For 
quasar B, we measure a H/3 FWHM of 2782±200 km s^^. 
The continuum luminosity at 5100A is 44.394 erg/s so 
logLBoi = 45.36. 

From IMcLure fc Dunlod (|2004[) . we adopt the black 
hole mass estimator 
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We thereby estimate black hole masses for quasar A 

and B such that log (^^%|^) = 8.46 and 7.80, respec- 

tiv ely. We estimate uncerta i nties of about 0.4 dex based 
on Vestcrgaard fc Petersonl (|2006|) . 

The Eddington luminosity, assuming a composition of 
pure hydrogen, is given by 
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where Mbh is the mass of the black hole, nip is the proton 
mass, and ctt is the Thompson scattering cross-section. 
Therefore we find the Eddington ratios I/Boi/.^Edd for 
quasars A and B to be 0.67 and 0.29 respectively. Com- 
paring to iShen et al.l ([200l for SDSS quasars in similar 
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ranges of redshift and LboI (see their Figure 12) , the Ed- 
dington ratio for quasar B is just 0.6 a above the mean 
in fogiBoi/^Edd, whereas for A the ratio is « 3(t high. 
This could be evidence that accretion rates are strongly 
boosted during cfose interactions among massive merging 
galaxies. 

4. MERGER SIMULATIONS 

We can further understand the properties of 
SDSS J1254+0846 via numerical simulations. Galaxy 
mergers may be a significant triggering mechanism for 
quasar activity, and there is growing interest in ver- 
ifying and understanding this connection more com- 
pletely. Since the majority of theoretical models (see, 
e.g., Wyithe & Loeb 2003, Volonteri et al. 2003, Hopkins 
et al. 2008) associate the most active phase of evolution 
(and thus most of the black hole growth) with nuclear 
coalescence, most quasars are expected to be hosted by 
systems where many of the telltale signs of interaction 
(disturbed morphology, tidal bridges and tails) no longer 
exist, or are difficult to find underneath the glare of the 
quasar. In man y cases, deep imaging at high spatial res- 
oluti on ([Dunlop et a l. 2003) but also high signal-to-noise 
(e.g.. lBennert et al]i2008l) is required to find evidence for 
these faint structures. 

In this context, SDSS J1254-I-0846, a pre-coalescence 
merger with two observed quasars, provides a unique 
opportunity to probe the early phases of the proposed 
merger/triggering mechanism. One of the most powerful 
insights into this system is via numerical simulation, i.e., 
designing numerical models of the current system which 
can be evolved, modified, and compared to the observed 
system. We have undertaken just such a task using nu- 
merical techniques that are extensively detailed in exist- 
ing literature (e.g., Springel et al. 2005, Cox et al. 2006, 
Hopkins et al. 2008, Jonsson et al. 2009). Briefly, initial 
equilibrium disk models are constructed to be represen- 
tative of disks at the appropriate redshift. These mod- 
els are then initialized on a prograde orbit and allowed 
to evolve using the N-body/SPH code Gadget (Springel 
2005) from a distant separation, through their interac- 
tion, to their eventual merger. The references provided 
above include extensive descriptions of these models, and 
their generic outcome. 

A representative result of such a numerical modeling 
experiment is shown in Figure [5l which displays a pro- 
grade parabolic galaxy merger with baryonic mass ratio 
2:1, viewed after the second close passage, but prior to 
the final coalescence. This model was selected owing to 
both its nuclear separation and the position and extent 
of the tidal features showing a remarkable resemblance 
to SDSS J1254+0846. On the other hand, the black hole 
masses are off by a factor of ^2 and their accretion rates 
by a factor of 5-10. 

Determining a suitable match proved to be a time- 
consuming endeavor which required the analysis of ^200 
merger simulations to isolated the orbits and orientations 
that best reproduce the observed tidal features, and the 
simulation of 8 additional mergers to perfect this match. 
In general, the symmetric tidal features place a relatively 
tight constraint on the relative orientation of the disk 
spins to that of the orbital plane. Specifically, the spin- 
orbit orientation is required to be less than ^^30 degrees. 
Furthermore, the relationship between the tidal features 



and the nuclear separation demanded a relatively large 
impact parameter (Rperi ^ '^Rdisc) such that it had a 
wide first passage, and a glancing second passage, prior 
to the final coalescence. Additional velocity information 
for both the nuclei and the tails would place even tighter 
constraints on these parameters. 

The observed tidal features offer a direct means to con- 
strain the orbital parameters. However, the set of ob- 
served galaxy properties - specifically the accreting black 
holes - provide insufficient information to uniquely de- 
termine the properties of the interacting galaxies. The 
matching experiment here only informs us that the pro- 
genitor spiral galaxies are required to be sufficiently large 
(i.e., scalelength of the stellar disc i?ciisc^4 kpc) , to pro- 
duce the length of the tidal tails, and that they must 
have contained pre-existing stellar bulges to match the 
black hole masses observed at this early merger stage. It 
might be possible to quantify the properties of the pro- 
genitors better with additional information about the ob- 
served system or using additional models, but we caution 
that the predicted black hole properties include assump- 
tions about the initial disk model, including the seed 
black hole masses, as well as black hole accretion that 
occurs well below our model resolution, and thus addi- 
tional model-matching is unwarranted until a larger sam- 
ple of observed galaxies exists. In particular, we should 
not expect (indeed, we should be skeptical of) a perfect 
match between the observed and model-predicted accre- 
tion rate onto the supermassive black holes (and thus 
also the luminosity of the quasars at any given time), 
because stochastic (unpredictable) accretion events ap- 
pear to turn ignite nuclear accretion activity at any 
given time along the merger sequence in ULIRGs (e.g., 
IVeil leux et al.ll20"09b[ ). The lower right panel of Figure [5] 
makes clear how noisy the accretion is expected to be. 

Given the success of our search for a model that 
matches many of the properties of the observed system, 
it is fair to ask what we learn from such an experiment. 
First, we have provided additional evidence that galaxy 
mergers are a plausible scenario for triggering quasar ac- 
tivity. In fact, nearly all quantities are extracted from 
the matched model. The orbit, the orientation, and the 
progenitor galaxies are fully consistent with, and even 
expected in, a merger-triggered scenario for quasar for- 
mation. Second, we have identified a case where the 
progenitor galaxies participating in the galaxy merger 
contain massive bulges and, hence, super-massive black 
holes. While most theoretical models do not currently 
make testable predictions about the abundance of such 
systems during the hierarchical growth of galaxies and 
black holes, a larger sample of such observed systems 
will certainly motivate additional investigation. 

5. A GRAVITATIONAL LENS? 

Given the identical redshift, similar colors, and the 
strong similarities in continuum slope and broad line pro- 
files between components A and B, is it possible that 
the pair is lensed? It seems highly unlikey for several 
reasons. First, the observed optical spectra are very dif- 
ferent, whereas gravitational lensing should be essentially 
achromatic. Second, the A/B flux ratio is unusually large 
for a lensed system. Third, to achieve the rather wide 
A/B separation, a massive lens is expected at interme- 
diate redshift, which should be sufficiently luminous as 
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to be evident in the images. Fourth, an intervening lens 
galaxy should produce some absorption signature in the 
quasar spectra. We discuss these four objections to a 
lens interpretation in turn. 

Spectral differences between quasar image components 
are common even in bona fide lenses (e.g., Wisotzki et 
al. 1993; Burud et al. 2002a,b; Oguri et al. 2005; Sluse 
et al. 2007). These differences are typically explained as 
the effect of either microlensing, or as light path time de- 
lays sampling intrinsic quasar spectral variability. Even 
with macrolensing only, anisotropy in the source may cre- 
ate spectroscopic differences along the slightly different 
sightlines (Green 2006, Perna & Keeton 2009). While 
particularly illuminating of source structure, such effects 
are expected to be much more subtle than those observed 
here. If the pair is indeed lensed, then microlensing re- 
mains the most likely explanation for the spectral differ- 
ences. However, microlensing-induced spectroscopic dif- 
ferences should primarily affect emission from the more 
compact emission regions of the source quasar - the con- 
tinuum, or perhaps the broad lines. The observed spec- 
troscopic differences are instead predominantly in the 
narrow lines, whose emission region is too large (hun- 
dreds to thou sands of parsecs; e.g., iMotta et all [200l 
iBennert et al.|[2002i ) to be affected by microlensing. 

Let us consider the second objection to a lens in- 
terpretation, the unusually large A/B flux ratio in 
SDSS J1254-h0846. While most known bona fide lens 
components indeed show smaller flux ratios, this is quite 
possibly a selection effect caused by flux limited surveys 
from which lens candidates are found. Furthermore, if 
the system were indeed lensed and microlensing were in- 
deed the cause of the spectroscopic differences described 
above, the unmicro- (but still macro-) lensed flux ratio is 
unlikely to be that observed in the broadband photome- 
try. Since microlensing is least likely to affect the narrow 
line region, the flux ratio of the narrow lines might more 
accurately represent the true macrolensed flux ratio. The 
A/B ratio of [O III] line flux above the continuum is only 
about 5.8±0.4. For comparison, the mean/median/mode 
of the A/B /-band (HST F814W) image flux ratio for 60 
lensed quasars in the CASTLE^l database (E. Falco, 
priv comm) is 5.4/2.6/1.5, but 8 of the 60 (13%) have 
ratios above 9. Therefore, the large observed image flux 
ratio in SDSS J1254+0846 does not on its own rule out 
a lens interpretation. 

The third objection is that a luminous lensing galaxy is 
expected to be visible. One possibility is that the lensing 
galaxy of this system happens to be the tidally disturbed 
system visible in our images. One can ask post-facto, 
how likely is it to find an interacting pair of massive 
spirals with spectacular tidal arms, well-centered near 
the mean position of the quasar images? If the spin of 
a galaxy is randomly aligned with the orbit, then only 
1 / 6th of the time is it aligned within 30 degrees of the 
orbital spin, and so only l/36th of the time would both 
be aligned, such as needed to produce the observed tails. 
It is difficult - perhaps fruitless - to attempt a further 
probability calculation for such an alignment, given the 
huge parent sample in which this exceedingly rare object 
was found, the complex selection effects, and the lack of 

CfA-Arizon Space Telescope Lens Survey information is 
hosted at http://www.cfa.harvard.edu/castles 



imaging and spectroscopy of comparable depths in large 
statistical samples. Clearly the most definitive test of 
the hypothesis of a tidally disturbed lens galaxy would 
be deep spectroscopy of the faint tails to determine the 
redshift of the associated stellar population. 

To further investigate the lens hypothesis, we have run 
lens models using GRAVLENS software (Keeton 2001)E3 
for a Singular Isothermal Sphere (henceforth SIS) at all 
intervening redshifts. For the observed total flux ratio of 
'--^9 we find reasonable velocity dispersions near 300 km/s 
for 0.15 < ziens < 0.3. From the /-band TuUy-Fisher re- 
lation (Masters et al. 2006), we expect a (Ty=300km/s 
galaxy to have absolute /-band magnitude —20.8 ± 0.2. 
For a redshift z ~ 0.22, the expected SDSS i band mag- 
nitude for such a galaxy is about 19.5, which is about the 
same as the fainter quasar. The expected position of the 
center of mass of the lens is determined in the model by 
the observed flux ratio. Larger A/B flux ratio means pro- 
portionally smaller distance from the lens to component 
A. In the Appendix, we have attempted to subtract the A 
and B nuclear point sources and determine the location, 
brightness and significance of any galaxy light between 
them. In summary, while there is some evidence for ex- 
tended emission around the quasar nuclei themselves, we 
find no evidence for significant extended emission with 
a centroid consistent with the expected lensing galaxy 
position. 

Fourth, a lensing galaxy might be expected to produce 
detectable signtaures in the quasar spectrum. We find in- 
sufficient contribution from stellar emission (either from 
a z = 0.44 host or from a putative intervening lensing 
galaxy) to create any detectable spectral features. For an 
absorber at a plausible z/ens = 0.22, the more commonly- 
detected optical/UV intervening absorption lines (such 
as Lya or the Mgll doublet near 2800A rest) are m 
the UV. Detection of the Call AA3934,3969A (H&K) 
or NalD AA5891,5897 absorption would be feasible, but 
since these are extremely rare (Wild & Hewett 2005), 
a non-detection here is not useful. Detailed inspection 
of the spectra of both A and B components yields no 
evidence for any significant absorption lines that might 
suggest an intervening (lensing) galaxy. 

X-rays sample intervening gas and dust in all 
phases. Our X-ray spectral fit including a ziens = 
0.22 absorption component, yields an upper limit of 
2.2x lO^^atoms/cm"^ from the A component spectrum 
for an unchanged best-fit continuum slope F = 2.01 ± 
0.08. The absorption upper limit from B is about three 
times weaker. 

On balance, we have considered several significant ob- 
jections to a lensing scenario, and we believe that de- 
spite the nearly identical redshifts, SDSS J1254-f 0846 A 
and B truly represent a binary quasar. Our simulations 
confirm what by all appearances is a merger with orbits 
fortuitously close to the plane of the sky, for which the 
very similar observed nuclear velocities are much more 
likely. 

6. CONCLUSIONS 

The quasars SDSS J125455.09-h084653.9 (A) and 
SDSS J125454.87-h084652.1 (B) are within 21kpc pro- 
jected transverse separation at their common redshift 



Software available at|http: //redfi ve.rutgers.edu/keeton/gravlens, 
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of z = 0.44, hosted by a galaxy merger showing clear 
tidal tail features. The quasar A/B flux ratio is nearly 
constant across all 5 SDSS bands, and they show a re- 
markably small ~2(T velocity difference of ~200 km s~^ . 
We find especially strong differences between their nar- 
row emission line equivalent widths, and their Edding- 
ton ratios. The spectroscopic features of A - in partic- 
ular weak [O III] [5007] with evidence for a blueshifted 
component - are associated with hig h accretion rates 
(jAoki et al.ll2005HKomossa et al.l200"a ) and outflows. We 
suggest that A is very strongly accreting. Given its weak, 
blueshifted [QUI], it would be a candidate for a broad 
absorption line quasar, verifiable with UV spectroscopy 
of the CIV region. A counterargument is that it has a 
normal X-ray/optica l flux ratio, which is rare i n BAL 
quasars (|Green et al.l [2"00I: GaU agher et al.ll2006[ ). 

The close coincidence of the positions, colors and red- 
shifts of the two quasar components raises the suspicion 
of lensing. We examine a variety of counterarguments, 
most prominently the strong optical spectroscopic emis- 
sion line differences, but also the large image flux ratio, 
and the absence of either emission or absorption signa- 
tures from an intervening lens galaxy. Deep spectroscopy 
of the tidal tails should prove interesting for a test of the 
lens model, but also for more detailed study of the stel- 
lar populations in this unique system. Although each 
counterargument to lensing may have known caveats, we 
find the overall evidence to be quite persuasive that the 
pair is indeed a binary quasar. Perhaps the strongest 
argument is simply the association of a binary quasar 
with a clear merger of two massive disk galaxies. While 
expected under the merger hypothesis for quasar trig- 
gering, we deem the coincidence of such a system with a 
lensing configuration to be exceedingly unlikely. 

Indeed, a simple explanation for the very similar nu- 
clear radial velocities, as suggested by the galaxy im- 
ages and numerical simulations, is that the merger is ap- 
parently along an orbit close to the plane of the sky. 
SDSS JI254-I-0846 may represent a rare system where 
the orientations of accretion disks in the quasar nuclei 
can be constrained by the system configuration as being 
close to our line-of-sight (modulo the unknown relative 
orientations between accretion disks and galaxy disks). 
Such an orientation would be consistent with the unob- 
scured. Type 1 spectroscopic nature of both quasars. 

There are strong advantages to studying spatially- 
resolved binary quasars (SRBQs) such as these. SRBQs 
can be particularly useful when chosen from well-defined 



parent samples. First, such samples probe ongoing merg- 
ers. Second, the spatial and velocity information, espe- 
cially when combined with well-resolved spectra provid- 
ing separate black hole mass estimates, offer more con- 
straints on the properties of the merging components and 
the evolution of the merger. We have found a good match 
via numerical merger simulations for the orbit, the ori- 
entation, and the galaxies in this system, showing that 
it is fully consistent with a merger-triggered scenario for 
quasar formation, where the progenitor galaxies already 
contain massive bulges. By selection, SRBQs are likely 
to be face-on, which makes them ideal for providing mor- 
phological constraints on merger models via e.g., follow- 
up with HST or ground-based adaptive optics. The use of 
uniform SRBQ parent samples further allows us to place 
these systems in their larger cosmological context, which 
is crucial if we are to understand the role of merger- 
triggered supermassive black hole accretion, and its re- 
lationship to galaxy evolution. 
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APPENDIX 
APPENDIX: SURFACE PHOTOMETRY 

We used the 2-dimensional galaxy- fitting program GALFIT (jPeng et al.l 12002^ to decompose the quasar and host- 
galaxy light of image components A and B and to detect and fit underlying extended features. GALFIT can si multa- 
neou s ly fit one or more obj ects in an image choosing from a library of functional forms (e.g., exponential, etc. iSersid 
119681 : IdeVaucouleursI 119481) . For convolution with the point-spread function (PSF) of the telescope optics, we first 
created PSF stars from stars within the image close to the location of the quasars. Most stars either had too few 
counts, especially in the PSF wings, or were saturated. We thus created an artificial Gaussian profile with a FWHM 
corresponding to that of non-saturated stars observed in the image. Qualitatively, the results are the same as using 
real star images, with the advantage of zero noise in the artificial PSF image. We find that >, 95% of flux is removed 
when we fit stars in the frame, with the residuals due primarily to ellipticity in the observed PSF. 

For both the r and the i-filters we subtracted the sky background and fitted two PSFs at the locations of quasar 
A and B. Simultaneously, we fitted and subtracted closely neighboring and bright objects. Due to the saturation of 
quasar A in our images, we fixed the fiux ratio of the PSFs of quasar A and B to the ratio given in SDSS-DR7 (using 
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"psiMag_r"; "psfMagJ"). We used the SDSS PSF magnitude for quasar B also to calibrate our results. Then, we fitted 
the ho st galaxies of quasar A an d quasar B a nd a potential underlying lensing galaxy with either a Ide VaucouleursI 
(|1948D profile or the more general ISersid ()1968[ ) profile 



E(r) = Scff exp 




(1) 



where Eeff is the pixel surface brightness at the effective rad i us rpa , and n is the Sersic index. In this generalized form, 
an exponential disk profile has n = 1, and a Ide Vaucouleurd (jl948l) profile has n = 4. In general, fitting a Sersic (1961) 
profile gives more flexibility to the fit, but also adds an additional free parameter to an already complicated fit, which 
can result in an unphysically large Sersic index or an unphysically small effective radius. For this reason, we decided 
to fix the Sersic index to either 4 or 1 (depending on the resulting x^) and set the minimum allowed effective radius to 
3 pixels (i.e., the minimum resolvable size given by the FWHM). We then chose the best fit based on the residuals and 
statistics. Note that we also constructed masks to exclude tidal structures during the fitting procedure; however, 
the effect of the latter on the results is negligible. 

In addition to fitting quasar A and quasar B with the PSF model, we followed these four approaches: (1) we fitted 
one "joint" host galaxy at a (starting) location in between quasar A and B; (2) we fitted two host galaxies at the 
locations of quasar A and B; (3) we fitted two host galaxies at the locations of quasar A and B plus another galaxy 
("lens") fixed at a position along the A - B line expected from a SIS lens mode^; (4) we fitted two host galaxies at 
the location of quasar A and B plus another galaxy ("lens") close to the location of quasar A. 

Only models 2 and 4 yield acceptable fits, with 4 slightly preferred by comparison of the reduced (i-e., taking into 
account the larger number of parameters in model 4) . Although it is easier to "hide" a lens galaxy near component A 
as in model 4, the lens model predicts that the brighter image ("A") is outside of the Einstein Radius and is further 
from the lens. The fainter image is closer to the lens and is interior to the Einstein Radius. Therefore, one result of 
our experiments is that host galaxies are required for both A and B. A single "joint" host is unacceptable: GALFIT 
instead preferred a host galaxy at the location of cither A or B. The second result is that a third extended component 
(a putative "lens" galaxy in addition to host galaxies for A and B) is pre ferred by the fits, but the best-fit position 
is not the one predicted by a simple SIS lens model. Fixing an additional Ide VaucouleursI (|1948D profile closer to the 
fainter quasar results in an unphysically tiny effective radius. Freeing the coordinates results in a second extended 
component (in addition to the host galaxy) at the location of quasar B, and one of the two components becomes either 
unphysically huge or vanishingly small. The third ("lens") galaxy is allowed only if it is much closer to quasar A. In 
Table [1] we summarize the results in both r and i filters for the best fit models, based on resi duals and s tatistics. 
Based on our ex perience with fitting quasar host galaxies (and simulations carried out; Bennert et al.|[2010t see also 
iKim et al1l2008l ). we conservatively estimate the uncertainties of the AGN luminosity to 0.2 mag and those of the host 
galaxies to 0.5 mag. 

Note that fitting this system is complicated and our results have to interpreted with caution. In general, the decom- 
position of complex images in multiple components is a difficult statistical challenge due the degeneracies involved, 
and the highly non-linear dependency of the likelihood on a large number of parameters. Decomposing quasar and 
host-galaxy light is already difficult. Here, the fitting is further complicated by the fact that we have two quasars close 
to each other, one of which is saturated, with possibly a merging host galaxy and/or tidal disturbances and/or another 
underlying galaxy. Keeping in mind these cautionary notes, we can conclude the following: There seem to be two host 
galaxies at the location of quasar A and B, not just a relaxed galaxy hosting two quasars. While there is no evidence 
for another galaxy close to the location of quasar B, we cannot exclude the presence of another galaxy close to the 
location of quasar A. However, based on the fiux ratios of QSO A and B, we would expect the lensing galaxy to fall 
closer to QSO B. Thus, our modelling rules out significant extended emission with a centroid that is consistent with 
the expected lensing galaxy position. Deep images of higher signal-to-noise and smaller PSF are warranted, either 
from HST or using adaptive optics from the ground. 
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Observed Wavelength (A) 

Fig. 1. — Optical spectra. The two components' spectra are plotted from 4500 - 9600 A, with the B component scaled up and shifted 
for clarity. Despite the factor of ~11 difference in flux normalization, the redshift, continuum and broad line shapes are all remarkably 
similar. The most striking difference is in the equivalent widths of the narrow emission lines; all the forbidden lines of [NeV], [Oil], [OIII] 
are relatively much stronger in B. Major emission line species are labeled along the top. Major atmostpheric absorption bands are marked 
with circumscribed crosses along the bottom. CCD gaps are evident in both spectra between ~6600 - 6700A. 



TABLE 1 
Results from GALFIT fitting 



Filter 

(1) 


PSF quasar A 

mag 


PSF quasar B 

mae; 


Host quasar A 

mae; 

U) 


Host quasar B 

mae; 

(•') 


"lens" 

mag 
(6) 


r 


17.15 


19.50 


17.52 


19.33 






17.15 


19.50 


17.98 


19.33 


18.23 


i 


17.07 


19.36 


18.44 


19.85 






17.07 


19.36 


18.75 


19.49 


18.37 



Note. — GALFIT host galaxy fit results. +Col. (1): Filter: SDSS r or SDSS i. Col. 
(2): SDSS PSF magnitude of quasar A. Col. (3): SDSS PSF magnitude of quasar B. Col. 
(4): Best-fit magnitude of quasar A host galaxy. Col. (5): Best-fit magnitude of quasar B 
host galaxy. Col. (6): Best-fit magnitude of third "lens" galaxy close to A. 
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Fig. 2. — Optical spectral features near H/3. In the upper panel, the two components' spectra are plotted from 6750 - 7600A, with the B 
component scaled up and shifted for clarity. The dashed curve is the (arbitrarily scaled) optical Fc II emission template of I Zw 1 (Boroson 
& Green 1992). Such iron multiplet emission appears to represent at best a minor contribution to either spectrum. The residual of 10.66*B 
— A is shown in the bottom panel. The fiat residuals highlight the similajrity of the redshifts, continuum shape, and broad line profiles. 
The apparent trough blueward of the stronger [O III] line is due to the broader, blueshifted profile of the line in A. 



SDSS J1254+0846: A Merging Binary Quasar 



13 




4SQ.0 WW eOQQ 



Observed. Wavelength .(A) 

Fic;. 3. — Optical spectral features near H7. In the upper panel, the two components' spectra are plotted from 4500 - 6500A, with the 
B component scaled up and shifted for clarity. The residual of 11.5*B - A is shown in the bottom panel. 
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Fig. 4. — Optical and X-ray images of SDSS J1254-I-0846. This optical image (1' on a side, TV up, E to the left) is the median of 
four 5min exposures obtained 2009 May 22 in l" seeing with an r band filter on the IMACS camera at the Magellan/ Baade telescope 
at Las Campanas Observatory in Chile. The two bright quasar nuclei are evident. The brighter A component (r = 17.5) was identified 
spectroscopically in the SDSS as a quasar. Discovery spectroscopy of the B component (r = 19.2) was obtained as part of our binary 
quasar survey. This follow-up IMACS image clearly reveals the tidal arms of a host galaxy merger. Inset: Our Chandra 16ksec ACIS-S 
X-ray image (same scale and orientation, but 0.5' on a side) shows the two nuclei, which both have typical fx/fopt and power-law spectral 
slopes (Green et al. 2009). There is no evidence for extended emission as might be expected from a host (or lensing-mass) group or cluster. 




Fig. 5. — Numerical simulation of a merger like SDSS J1254+0846. The left panels show the gas distribution, and the middle panels the stars in the xy and xz planes, for this model 
of a prograde merger of two massive disks. The right panels show the model star formation rates, nuclear black hole masses and accretion rates for the two components. The epoch 
displayed in the images at left, 2.3 Gyr, is marked with red diamondsin the plots to the right, at which the nuclear separation, the position and extent of the tidal features (in the xy 
plane), the black hole masses, and their accretion rates resemble those of SDSS J1254+0846. 



